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ABSTRACT 

We report the discovery of 854 ultra diffuse galaxies (UDGs) in the Coma cluster using deep R band 
images, with partial R, i, and Ha band coverage, obtained with the Subaru telescope. Many of them 
(332) are Milky Way-sized with ver y large effective rad i i of > 1.5 kpc. This study was motivated by 
the recent discovery of 47 UDGs by van Dokkum et al. (2015a); our discovery suggests > 1,000 UDGs 
after accounting for the smaller Subaru field (4.1 degree^; about 1/2 of Dragonfly). The new Subaru 
UDGs show a distribution concentrated around the cluster center, strongly suggesting that the great 
majority are (likely longtime) cluster members. They are a passively evolving population, lying along 
the red sequence in the color-magnitude diagram with no signature of Ha emission. Star formation 
was, therefore, quenched in the past. They have exponential light profiles, effective radii ~ 800 pc- 
5kpc, effective surface brightnesses ^J,e{R) =25-28 magarcsec”^, and stellar masses ~ 1 x IO^Mq - 
5 X 1O®M0. There is also a population of nucleated UDGs. Some MW-sized UDGs appear closer to 
the cluster center than previously reported; their survival in the strong tidal field, despite their large 
sizes, possibly indicates a large dark matter fraction protecting the diffuse stellar component. The 
indicated baryon fraction < 1% is less than the cosmic average, and thus the gas must have been 
removed (from the possibly massive dark halo). The UDG population is elevated in the Coma cluster 
compared to the field, indicating that the gas removal mechanism is related primarily to the cluster 
environment. 

Subject headings: galaxies: clusters: individual (Coma) - galaxies: evolution - galaxies: structure 


1. INTRODUCTION 


This study is motivated by the discovery of 47 ul¬ 
tra diffuse galaxies (U DGs) in the Coma cluster by 


van 


Dokkum et al. ([201ba | using the Drag onfly Telescope Ar¬ 


ray (| Abraham &: van Dokkum 2014 hereafter Dragon¬ 
fly) . This unexpected discovery revealed a new popu- 
latoin of low surface brightness (SB) galaxies. Indeed, 
their central SBs are very low 24-26 magarcsec”^ in g- 
band and their median stellar mass is only ~ 6 x IO^Mq, 
despite their effective radii re = 1.5-4.6 kpc being as large 
as those of L* galaxies ( e.g., ^ 3.6 kpc for t he Milky Way 
(MW), calcul ated from |Rix & Bovy|2013|. |van Dokkum 

UDGs probably have 
gh dark matter fractions as they have survived in 


the strong tidal field of the cluster. 

Dragonfly is optimized to efficiently discover faint 
structures over a large field of view, but has relatively 
poor spatial resolution with seeing and pixel scales of 
^ 6" and 2.8", respectively. The above properties of 
the Dragonfly UDGs were derived after their discovery 
from archival Canada France Hawaii Telescope imaging. 
Follow-up studies are needed to understand their nature 
and origin, as well as their relationsh ip to the cluster en¬ 
vironment (Boselli &: Gavazzi 2014 for review) and to 
other more studied galaxy populations, such as dwarfs 
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and low SBs, in clu sters (e.g., Binggeli fc Cameronj 1991 [ 


119911 Ulmer et al. 119961 2L)11| Adami et al. 

2006b 2009[ |Lieder et al.||2012 ]Ferrarese et al.||2012 ). 

Optical telescopes of larger aperture a re advantageous 
for a r esolved study of this population. [Yamanoi et al.| 
(2012) used the Subaru Pri me Focus Camera (Suprime- 
Cam; Miyazaki et al.|2002 1 on the Subaru telescope and 
derived a galaxy luminosity function down to Mr ^ —10 
in Coma. Their three 34' x 27' fields include nine Drag¬ 
onfly UDGs. All of the nine were easily found in their 
catalog, being resolved spatially in the images. There¬ 
fore, Subaru imaging can identify this new galaxy pop¬ 
ulation efficiently and permits an investigation of their 
internal properties. Several ar chival Subaru images are 
available for the Coma clust er (|Yagi et ah 2007 Yoshida 


Bothun et al. 


et al.|2008 Yagi et al.|2010 (Jkabe et al.|2010'2014|). In 

this Letter, we use the archival Subaru data and report 

the discovery of 854 UDGs, implying ~ 1000 UDGs after 
scaling for the Dragonfly field-of- view. __ 

[20dd| 

:h cor- 


We adopted m — M = 35.05 (Kavelaars et al. 
as the distance modulus of the Cdma cluster, whic 
responds to an angular diameter distance of 97.5 Mpc 
(1" = 0.47 kpc) R The full catalog of the Subaru UDGs 
will be publisheu in Yagi et al. (in preparation). We use 
the AB-magnitude system in this work. 

2. DATA 

The raw R band images from the Supri me-Cam were 
obtai ned from the Subaru data archive (Baba et al.j 
2002). Suprime-Cam has a mosaic of ten 2048 x 4096 
GGDs and covers a wide field of 34' x 27' with a pixel scale 
of 0.202". An eighteen-pointing mosaic with Suprime- 


® We adopted the C osmological param eters of {ho,Q,M, 
O;,) = (0.71, 0.27,0.73) from [Larson et al.|||20ll|. 
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Fig. 1. — T he 2.86 deg x2.90deg 4.87 x 4.94Mpc^) area centered on the Coma cluster, the same area as in Figure 1 of I van Dokkum 

H ( [2015a| |. 'a) Image from the Digitized Sky Survey. The white borders show the 18 fields covered in the Subaru R band_^)tob^^^aU 
, which have the total area of 4.1 degre e^, about 1/2 of t he Dragonfly coverage. Red indicates the area analyzed by |Yamano ^^^ 
Yellow outlines the area analyzed by l^^giet al. l|20Tol using the Subaru B, R, Ha, i bands. Cyan indic ates the area m figure 
rne center of the cluster (aj2ooOi'5j200o)=(12:59:4:^.8,+zY:t)»:l4) is marked with a green cross ||White et al.|1993||. (b) The same area as in 
(a), showing the distribution of the 854 Subaru UDGs (circles). The MW-sized UDGs, with large effective radii (> 1.5 kpc), are shown in 
blue. The Subaru field coverage in R is enclosed with the solid line. The 47 Dragonfiy UDGs are indicated with red crosses. 


Okabe et al. (2014) and covered about 


Cam was taken b }JU 

4.1 deg^ (Figure [l|. The seeing was 0.6-0.8", typically 
0.7". Integration times for the 18 fields were not the 
same, resulting in variations in background noise, i.e., 
28.3-28.7 mag arcsec“^ (Icr) in a 2" aperture (equivalent 
to 30.0-30.4 mag arcsec~^ in a 10" ape rture, ^ 1 mag 


deeper than van Dokkum et al. 2015a). The very cen¬ 
tral field has a higher variation of 27.8 mag arcsec“^ since 
the field is contaminated by the outer envelope of bright 
galaxies. 


The data were reduced in a standard way (Yagi et al. 
2002 2010). We used self-sky flat images, subtracted sky 


background locally in each small gri d (256 x 256 pixel^ 


51.7" X 51.7"), used the WCSTools (Mink 


2002) for as¬ 


trometry, a nd applied a photometric calibration (Yagi 
et al.||2013|) usin g the Sloan Digit al Sky Survey (SDSS)- 
111 DKO catalog (Ahn et al. 12012 ). The grid size for the 


background subtraction was larger than the expected size 
of UDGs (< 30" ~ 15kpc). The Galactic extinction in 
R band vari es from 0.016 to 0.031 mag across the 18- 
field mosaic (Schlafly & Finkbeiner 2011). We adopted 
the Galactic extinction value at the center of each field 
and neglected variation within each field. The final pho¬ 
tometric error is < 0.1 mag. More details of the data 
reduction procedure will be presented in Yagi et al. (in 
preparation). In addition, we used the Supri me-Cam B, 
i, and Hct reduced images (see F igure fT]) by |Yagi et al. 
(2010) and Yamanoi et al. (2012‘ 


We also analyzed a control held for comparison. The 
R band data of one Suprime-Cam pointing, 1/18 of the 
Coma field, were taken from the Subaru Deep Field 


(SDF) project (Kashikawa et al. 2004). We used only 
a part of the raw SDF exposures taker! in June 2008 to 
make the background noise comparable to that in the 
Coma fields. The Icr background noise in a 2" aperture 
is 28.6 mag arcsec“^. For consistency we started from the 
raw data and matched data reduction parameters. 

3. IDENTIFICATION 

Our goal is to find UDGs in the Subaru images. Forty 
of the 47 UDGs discovered by Dragonfly are within the 
Subaru R band covera ge based on their coordinates (van 
Dokkum et al. 2015a). All were detected significantly 
(with the faintest one, DF27, off by 12.5" from the pub¬ 
lished coordinate) and their structures were resolved in 
the Subaru images. The detection threshold was approx¬ 
imately 27.3 magarcsec”^ in the R band. We describe 
our selection procedure for the final catalog of 854 UDGs 
in the Coma cluster. We found no counterparts in the 
control field. 


We ran SExtractor (version 2.19.5; Bertin & Arnouts 
1996) on individual mosaic frames with a hxed de- 
tection threshold of 27.3 magarcsec”^ in R. We re¬ 
moved a hrst set of spurious detections using SExtrac- 
tor’s ’’FLAGS < 4” and ’’PETRO.RADIUS > 0”. This 
initial catalog had 2,627,495 objects, including dupli¬ 
cates in the overlap regions of adjacent mosaic frames 
(^ 30 %). We used the Dragonfly UDGs as the fidu¬ 
cial set in adjusting parameters for selection of UDG 
candidates, but could not use exactly the same selec¬ 
tion criteria as van Dokkum et al. (|2015al due to the 
difference in image quality. We applied constraints on 
R magnitude and size, ”18 < MAG_AUTO < 26” and 
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Fig. 2.— Subaru BRi color image of the ~ 6' X 6' region (~ 170 X 170kpc^ region at d = 97.7Mpc), shown in cyan in Figure [l^. The 
Dragonfly and Subaru UDGs are marked respectively with yellow and green circles with a diameter of 20” (~ 9.5 kpc). 


”FWHM(Gaussian) > 4"” (i.e., all Dragonfly UDGs sat¬ 
isfy this condition), which left 7,362 objects. 

The reported effective radius of th e Dragonfly UDG s 
is re > 3.2" (using from GALFIT; Peng et al.||2002). 
However, in the resolved Subaru images we found that an 
alternative constraint, SExtractor’s > 1.5", captures 
all the Dragonfly UDGs. We therefore used re > 1.5" and 
a mean SB of (^(re)) > 24 to choose UDG candidates. 
[Note that we found that from SExtractor and GAL¬ 
FIT were occasionally very inconsistent; we use SExtrac¬ 
tor’s re for identification and GALFIT’s re for studies 
of structural properties.] We excluded objects with high 
central concentrations (mostly high-z galaxies) by remov¬ 
ing those whose mean SB within re deviates largely from 
the SB at re. This constraint, Re — (/i(Fe)) < 0.8, left 
1,779 candidates. 


The final step was removal of spurious objects by vi¬ 
sual inspection. Most spurious detections were due to 
the crowding in the cluster, such as faint tidal tails and 
galaxy blending, as well as distant edge-on disk galaxies, 
artifacts at image edges, and optical ghosts. To minimize 
human error, the four authors separately went through 
all postage stamp images. After this step and removal of 
duplications based on their coordinates, 854 UDG candi¬ 
dates were left on which at least three of us agreed. The 
full catalog will be published by Yagi et al. (in prepara¬ 
tion). 


4. ULTRA-DIFFUSE GALAXY CANDIDATES 

The 854 UDGs candidates from Subaru are visually 
comparable to the Dragonfly UDGs. Figure shows a 
sample 6' x 6' field, showing the Subaru (green circles) 
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and Dragonfly UDGs (yellow). Their low SBs are evident 
compared to the surrounding galaxies, including ma¬ 
jor galaxies in the cluster and distant background ones. 
Their large sizes are also clear when compared to the 
20" diameter of the circles (~ 9.5 kpc at d = 97.5 Mpc). 
The greater number of detections, compared to Dragon¬ 
fly, may be due to the superior seeing (less blending) and 
higher signal-to-noise ratio. 

The majority of the 854 candidates are most likely 
UDGs in the Goma cluster. One of them has been 


/io(R) ranges arou nd 23-26 mag arcsec Th ese results 


spectroscopically conf irmed as a cluster member (van 
Dokkum et al. 2015bI. The control SDF field has vir- 
tually no counterparts - only thirteen were left after the 
SExtractor-based selection, twelve of which were obvi¬ 
ous image artifacts or tails of bright galaxies. The last 
one appeared to be a blend of multiple objects. Hence, 
contamination by non-cluster members is rare and neg¬ 
ligible. Note, however, that some negligible number of 
contaminations might still exist. For example, the third 
object from the top in Figure may be a background 
spiral galaxy. Despite this significantly increased sam¬ 
ple, the UDGs are still a min or population in the Coma 
cluster (|Yamanoi et al.||2012|). 

In the literature, we found that many of the Subaru 
UDG s had been cataloged albeit as more compact ob¬ 
jects; Adami et al. (2006a ) found 248 of 309 t hat lay 
within their coverage, and Yamanoi et al. (2012) 232 of 
240. A mong them, only 17 were classified as low SB 


galaxies (Adami et al. 
low SBs were revealed 


2006b). Their large extents and 
or the first time in this study. We 


note that eleven out of the 12 D ragonfly UDGs within 


their field were also cataloged in Adami et al.| (2006a ) 
but none were classified as low SB ([Adami et al. 2006b). 


5. STRUCTURAL PARAMETERS 

The GALFIT package was used to measure the struc¬ 
tural parameters of the Subaru UDGs. The fits were 
made with a single Sersic profile ( Sersic||1968 1 with sky 

l-\ r\ ^ GT \A/r^ t- i <-1 w L<’it**!-•»» ' n ^4~ t-\ 4~-rt 


background fixed. We used SExtractor s segmentation 
images to mask surrounding objects. Seventy-nine of 
the 854 required additional manual masks to exclude a 
bright compact object(s) within their boundaries - inter¬ 
estingly, sixty-seven of these appear to have compact nu¬ 
clei at their centers. We judged the fits acceptable based 
on the goodness-of-fit {xt < Ij 75 objects were thus ex¬ 
cluded) and consistency between GALFIT and SExtrac¬ 
tor measurements {r^ from GALFIT and SExtractor con¬ 
sistent within a factor of 3; 11 removed). In this section, 
we use the sample of 768 objects with good GALFIT re¬ 
sults, out of which 332 have G ALFIT’s > 1.5 kpc (i.e .. 


MW-sized UDGs as defined in van Dokkum et al. 2015a). 
We refer to the former full set of galaxies as UDGs, and 
the latter as MW-sized UDGs. Figure!^ shows some ex¬ 
amples of GALFIT results for the UDGs of lowest SB, of 
largest-size, and with a compact nucleus. 

Figure l^a), (b), and (c) show histograms of Sersic in¬ 
dex (n), axis ratio (b/a), and central SB (^o(R))- For 


are consistent with van Dokkum et al. (2015a) when the 
difference in the adopted bands, SDSs g and Subaru R, 
is taken into account (roughly g — R ~0.8 mag for the 
red-sequence in Coma). 

The Subaru UDGs are likely an extension of normal 
and dwarf galaxy populations and are not, on their own, 
a distinct population. Figure |^d) shows the properties 
of the UDGs (crosses) w ith respect to normal g alaxies 
in Coma (circles; from |Komiyama et al. 2002). The 
apparent R magnitude of the UDGs is 18 to 24 liiag, in¬ 
dicating an absolute magnitude of about -12 to -16 mag 
at the Coma distance. The smallest and faintest UDGs 
(e.g., Te < Ikpc and Mr < —12 in Figure]^) over¬ 
lap with the largest and brightest dwarf galaxies and 
share some properties in common with them ( e.g., the 
exponential profile, nucleated population; see [Tolsto; ' 




et al.|2009 McConnachie||20 12 jBoselli & Gavazzip014p. 

Dotted lines represent constant SBs (^eS) from 23 to 29 

magarcsec“^ with a 1 magarcsec“^ interval, assuming 
an exponential profile. The average SB of the Subaru 
UDGs is distributed from about 25 mag arcsec”^ (i.e., 
a cut-off due to the selection) to 28 mag arcsec”^ (due 
to the detection limit; this lower boundary is lower than 
the pix-to-pix detection limit, because the UDGs are ex¬ 
tended). 

The absolute magnitudes correspond to stellar masses 
of 1 X lO^M0-5 X 10*Mq if we adopt a mass-to-light ratio 
of M/Lr ~ 3. Note the M/Lr varies by a factor of ~ 2 
for ages of 4-12 Gyr and metallicities betwe en 0.2-1.0 
solar based on calculations using Starburst99 (Leitherer 
et aHl999 ), a single starburst, and a Kroupa initial mass 
function. 

6. A PASSIVELY-EVOLVING POPULATION 

The UDGs are distributed widely over the entire area 
of the cluster with a concentration toward its center (Fig¬ 
ure This spatial correlation also supports the as¬ 

sumption that the great majority are cluster members. 
Figure nearly covers the virial radius o f the cluster 
(^ 2.8 Mpc; ~ 1.7deg; Kubo et al. 2007), and reveals 
their relatively symmetric distribution around the cen¬ 
ter with a potential elongation toward the south west 
(roughly toward NGC 4839). This symmetric, wide¬ 
spread distribution may indicate their long history within 
the cluster. 

The UDGs closely follow the red-sequence of a 
passively-evolving galaxy population on the color- 
magnitude diagram . 232 UDGs are in the catalog of |Ya-| 
maimi et al. (2012) with both B and R photometry. Fig- 
ureshows their distribution (green). The comparison 
data (red and blue) show other cluster member gala xies, 
as well as backgrou nd galaxies, and are also from [YaT] 


manoi et al.| (2012). T he red-sequence is evi dent, and 


the solid line is a fit by Yamanoi et al. (2012). Clearly, 


the UDGs lie along this red-sequence, and their B-R 
colors are around 0.8-1.0 mag. This is similar to the 


sic indices (n) = 0.9-1.0 indicate an exponential pro- 

ters 1 

Adami et al.||2009 

Ulmer et al. 

2011 

Lieder et al. 

file. The distributions of axis ratio, as well as its av- 

2012). 


erage {b/a) = 0.7-0.8, are skewed toward a large value; 
therefore, this UDG sample does not consist of randomly- 
oriented thin-disk galaxies in a statistical sense (which 
would skew their distribution toward a low b/a). The 


"No significant Ha excess was found in UDGs. 217 
UDGs are within the Subaru Ha coverage (yellow in 
Figure [^), which was designed to det ect faint Ha emis - 
sion around Coma member galaxies (Yagi et al. [2010 ). 
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Largest Size Lowest Surface Brightness Nucleated 



Fig. 3. — Examples of GALFIT results drawn from the groups of largest-size UDGs, lowest surface-brightness UDGs, and nucleated 
UDGs. 



[mag/arcsec'^ ] 



Fig. 4.— Structural properties of UDGs. (a) Histograms of Sersic index n, (b) axis ratio 6/a, and (c) central SB with their 

medians, averages, and standard deviations. Black lines are for all 854 UDGs, while blue are for 332 MW-sized UDGs alone, (d) Effective 
radius vs. R magnitude. The parameters of the UDGs (crosses; red for the D ragonfly UDGs) are derived with GALFIT . Normal galaxies 
(circl es) -spestroscopically-confirmed Coma members ( [Mobasher et ar]|2001| | — are also plotted for comparison (from [Komiyama et ah] 
|2002[ with the conversion R(AB)-R(Vega)=0.21). Dotted, diagonal lines show constant SBs (/^es) from 23 to 29 mag arcsec ^ with a 1 
mag arcsec~^ interval for the case of an exponential profile (note /xq = Me — 1-82 for n = 1). The gap between the normal galaxies and 
UDGs is due to selection effects. Horizontal lines show rg of PSF with a FWHM of 1.5 arcsec (Komiyama et al. 2002) and a FWHM of 
0.7 arcsec (this study). 
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Therefore, the UDGs are not forming stars at the cur¬ 
rent epoch, as expected for passively-evolving galaxies. 


7. DISCUSSIONS 

We report the discovery of ~ 1000 UDGs in the Coma 
cluster, about 40% of which are MW-sized. The new 
UDG sample is by no means complete, but alr eady 
contains 10-20 times more than previously known (van 


Dokkum et al. 2015a). None of the UDGs show a sig- 
nature of tidal distortions; this is our selection criteria, 
but indicates that this sample of UDGs are not likely 
recently-disrupted tidal debris. 

The UDG population, compared to brighter galaxies, 
is elevated in the Coma cluster, although a small num¬ 
ber of large, low-SB galaxies are known in the field (jPal- 
canton et al.||1997l |Impey et al.||20011 1 Burkholder et al. 
200ip . t the UP G-to-brighter galaxy number ratio in 
Coma were common in the field, the expected UDG pop¬ 
ulation would be implausibly large, > 10® within 100 
Mpc of the MW. To obtain this rough estimate, we used , 
as a reference, galaxies in the SPSS (Ahn et al. 2012) 
within 16 < r < 17 mag (—19 < Mr < Ifl in ab¬ 
solute magnitude) at t he cluster’s redshift of 0.013-0.033 
(Mobasher et al. 20011. The number of reference galaxies 
in the held was e s tiiriat ed from the luminosity function 
of Blanton et al. (2001). The estimated number in the 
field, > 10®, is crude, but seems too large compared to 
the small number discovered so far. The cluster environ- 
ment must play a role in thei r formation and evolution. 


van Dokkum et al. (2015a) speculated that the MW- 
sized UDGs might be a dark-matter (DM) dominated 
population in order for them to survive in the strong 
tidal held around the cluster core. In fact, the Dragonhy 
UDGs spatially avoided the central r ~300 pc region as 
if the ones there had been tidally disrupted; this appar¬ 
ent disruption w as used to constrain the D M fraction (as 
large as > 98%; van Dokkum et al.|2015a ). Surprisingly, 
we found UDGs even closer to the core (Figure]^), and 
the closest one is MW-sized only about 3' (^ 85kpc) 


away on the sky. Eleven UDGs were found within a ra¬ 
dius of 5' (^ 141 kpc). These detections were, of course, 
not complete due to the high background emission there, 
and their apparent proximity may result from a chance 
coincidence along a line of sight. If any of them are 
within ^ 100-150 kpc from the core, an even larger D M 
mass than the estimate by van Dokkum et al. (2015a|) is 
necessary, and the baryon traction within a tidal radius 
should be < 1%. This is below the cosmic average, and 
therefore, the baryons must have been removed from the 
possibly very deep DM potential. 

The possible removal of the gas, and quenching of 
star formation (SF), are consistent with their red color 
and clustering in Coma (indicating their longevity within 
the cluster). The red-sequence can be produced by a 
metallicity-sequence if galaxi es have been evolving pas - 
sively since SF was quenched (Kodama & Arimoto|1997 ). 
Physical pr ocesses often suggested tor the queiich in- 
clude (see Boselli & Gavazzi 2014, for review): (a) 
blow out of gas d ue to galactic winds from supernovae 
or AG N activities (Dekel & Silk||19^ Arimoto fc Yosh iil 
T987| ), (b) ram-pressure strippmg (|Gum^Gott|jl9^ , 


(c) tidal-interaction and harassmenrlIVloore et ai.|199bh 
and (d) st arvation due to the cessation of gas intall ([ban 


son et al. 1980). The elevated population in the cluster 
indicates that environmentally-driven mechanisms, such 
as (a), (b), and (c) are the most likely solutions [(a) may 
occur if SF is induced by (b) or (c).]. 
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Fig. 5.— Color-magnitude diagram using B and R band photometry. The green points are 232 UDGs observed both i n B and R with 
Subar u (the Dragonfly UDGs are circled), and the red and blue are red and blue galaxies taken from the Comal field of|Yamanoi et al.| 
1(20121 which includes cluster members as well as background galaxies. Due to saturation, most giant galaxies are not included, but the 
red-sequence is evident. The UDGs clearly follow the red-sequence population of the Goma cluster. 
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